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ABSTRACT 


The motion of a large space structure, such as a space station, is described by a large 
number of coupled, second order differential equations. To effectively control this 
structure, a mathematical model is required. Both the mathematical model developed 
directly from the physics of the structure, and the simplified model developed with modal 
analysis are of extremely high dimension. A reduced order model is therefore required 
in order to design a control system for the structure. 

A straightforward approach to the control problem is taken by using linear quad- 
ratic optimal control techniques to determine the reduced order control solution for the 
truncated modal model. The effects of reduced order modeling on the control of the 
space station will be evaluated by observing the response of the closed loop system to 
several disturbances. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. BACKGROUND 
A large space structure, such as a space station, presents numerous control problems 


to the engineer. One such problem involves the control of unwanted vibrations. These 
vibrations must be controlled to prevent the disruption of delicate scientific experiments 
and maintain the stability of the structure. 

The lightweight materials used 1n the construction of the space station and its large 
size combine to form a flexible, ightly damped structure that will vibrate for a consid- 
erable length of time when disturbed. This thesis addresses the problem of actively 


controlling these vibrations. A representation of a dual keel space station, courtesy of 


McDonnell Douglas Astronautics, is presented in Figure 1. 
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Figure 1. Representation of a Dual Keel Space Station. 


B. PROBLEM STATEMENT 

To solve the vibration control problem, several steps are required. First, a math- 
ematical model must be developed that describes the behavior of the system over time. 
Modal analysis develops the modal model representing the system as a set of uncoupled, 
simultaneous, second order differential equations. The order of this model is high, 
making control design and implementation difficult. A reduced order model can be 
easily generated by truncating the modal equations. 

Second, a control system must be designed using the mathematical model. The 
control system design is performed using optimal control techniques, that is, a control 
solution is found that will minimize a performance or cost function of the structure. In 
this case, the cost function will be the vibrational energy of the structure plus a term 
representing the control energy. This cost is then minimized through the development 
of an optimal feedback gain matrix. 

Third, the controlled model must be simulated to determine the system response to 
a disturbance applied to various points on the structure. This simulation will utilize 
various optimal gain matrices based on a range of reduced order models. The effects of 
reduced order modeling on the control of the space station will be evaluated by observ- 
ing the response of the closed loop system to several disturbances. 

Finallv, conclusions will be presented based on the results, and recommendations 


will be made for further areas of research. 


C. ORGANIZATION 

In Chapter II, the model of a space station will be developed. The modal model is 
developed and discretized to yield the discrete-time state equations used in the simu- 
lation. The data for this model was provided by the McDonnell Douglas Astronautics 
Company. 

The desired performance function will be obtained in Chapter III, along with the 
optimal gain matrix necessary to minimize this function. Together, these provide the 
basis for the reduced order control solution. 

Chapter IV presents the computer simulation of the model and the reduced order 
control solution. The results of reduced order control on the system will be evaluated. 

Conclusions based on the simulation results are presented in Chapter V, as well as 
recommendations for areas of future investigation. 


ll. THE MATHEMATICAL MODEL 


A. INTRODUCTION 

A mathematical model describing the motion of the space station is required for ef- 
fective control. Because of its size and construction, the space station can be considered 
to be a lightly damped, vibrating structure which, when disturbed by an external force, 
may vibrate for a considerable length of time. 

The space station can be modeled as a finite number of discrete masses connected 
by springs and dashpots. This is a complex description of the system’s motion and 1s 
quite difficult to work with. By expressing the equations of motion in terms of the 
structure's natural modes of vibration a simpler model will be obtained. 

The natural modes of the system will form the basis for the continuous-time math- 
ematical model, developed in section B, expressed in terms of an uncoupled set of second 
order differential equations. The discrete-time model will be derived in section C from 


the results of section B. 


B. THE MODAL MODEL 

The space station is a lightly damped structure consisting of numerous natural 
modes of vibration. The structure can be modeled as a system of discrete masses con- 
nected by springs and dashpots. In mechanical svstems, the springs represent the 
suffness factor and the dashpots represent the damping factor of the system. The dis- 
placement of the masses can be described by a second order matrix differential equation 


of motion: 
M &(2) + on K a(1) + K a(n) = Fl) (2.1) 


where: 
e q is the generalized coordinate vector 


e M 1s the diagonal system mass matrix 


d Pe a 
‘ ion K is the structural damping term 


e dis the damping constant 
® q, is the frequency of oscillation of the system 


e K is the symmetric system stiffness matrix 


e F(t) is the system forcing function. 


This equation represents a system of simultaneous, second order differential equations 
that are coupled through the K matrix. 

Equation (2.1) can be decoupled by expressing q in terms of the natural modes of 
vibration. This is known as the process of modal analysis. The system is then repres- 
ented by a set of independent differential equations that can be treated individually. 
Hurty and Rubinstein (Ref. 1] and Meirovitch (Ref. 2] outline the modal approach that 
forms the basis for the model developed in this section. 

The first step in the model’s development is the solution of the undamped, homo- 
geneous form of Equation (2.1): 


M (1) + K q(t) =0. (2.2) 


The solution to this equation can be found in any elementary differential equations 
textbook and can be written as: 


q(t) = Ax sin(wr + ©) (2,3) 
and 
q(t) = —Axw’ sin(wr + ©). (2.4) 
Substituting these expressions into Equation (2.1) and solving: 
{[K — o M}x=0 (2.5) 
or 
Kx = w*Mx. (2.6) 


Equation (2.5) is an eigenequation with m combinations of x and w as solutions. 
Grouping the individual solutions: 


X = [x,x% .. XJ) (2.7) 
Q? = diag| w3),059, -- an |, (2.8) 
all solutions can be found by solving the matrix eigenvalue problem: 


KX = Q’MX (2.9) 


where {2? is the system natural frequency matrix and X is the system eigenvector or 
modal matrix. The individual elements of 2 are referred to as the natural frequencies, 
or eigenvalues of the system, and the columns of X are referred to as the natural mode 
shapes or the eigenvectors. Several properties of the eigenvector matrix are useful and 
can be developed by premultiplying Equation (2.9) by Xt: 


X'KX = Q°?X'MX. (2.10) 
The eigenvectors can be normalized: 
X'MX =I (2.11) 
where | is the identity matrix. From Equation (2.10), it follows: 
X'KX = 9°71 = Q’ (2312) 


where 92? is a diagonal matrix. 
The equations of motion can be uncoupled through a linear transformation of the 
coordinate system [Ref. 3]: 


A 


qi) = > xin) = Xn(0 (2.13) 


i= 


where q(t) and 7(t) are two different sets of generalized coordinate systems and 
e X 1s the modal matrix 
¢ nis the maximum number of degrees of freedom 


e® n(t) is the transformed coordinate vector or the modal amplitude vector. 


Applying the transformation to the system results in: 
M Xn(t) + 7 KXn(t) + K Xn(t) = F(d). (2.14) 
Multiplying both sides of Equation (2.14) by XT: 
X™M Xi(1) + = XTKXn(t) + X™K Xn(t) = X™F(0). (2.15) 
Applying Equations (2.11) and (2.12), noting that for wideband excitation, w,% o,; 


n+dQn+Q’y = X'F (2.16) 


where {2 has been established as a diagonal matrix. Equation (2.16) is the modal model, 
representing an uncoupled set of second order differential equations describing the mo- 


tion of the structure in terms of its natural modes of vibration. 


C. DISCRETE-TIME MODEL 
The discrete-time, state space model is found by solving the continuous-time 
equations obtained in Section B. Solving the modal system of equations for the i sol- 
ution determines the solution for all the individual second order differential equations 
of motion. The solution to this single equation will be used in the computer simulation. 
The :* equation of motion 1s expressed as: 


ni!) a dwoin{t) Ir wand!) = x} F(t) (2.17) 


where: 
® xi is the transpose of the i* mode shape vector 


e F(t) is the torquing force applied at a point. 


The homogeneous solution of this second order differential equation is: 
ge Pre 
n{t)=C,e ° 2 ' cos(ut)+C,e ~~ 2 * sin(ue) (2.18) 
where 


8 / Aw?, oa d°w?, 


= 5 ; (2.19) 


By defining 


d : 
y= —2 (2.20) 


and 


2 =o) =e (2.21) 


then substituting Equations (2.20) and (2.21) into Equation (2.18), the solution can be 
written: 


n{t}=C,e” cos(wgt)+ Ce” sin(wgt) (2.22) 


and 


nit) = (CO Ciy)je™ cos(wat) —(Cywg + Coy)e™ sin(we?). (2.23) 
Att = 0, Equations (2.22) and (2.23) become: 
n(0) = C, (2.24) 
and 
n(0) = C,wg— Cy. (2.25) 


Solving for C, and C;;: 


C ae Ba 
= ! a 
ce a - n(0) aii 


The homogeneous solution to Equation (2.17) can then be written in terms of C, and 
eC: 


mo) | eT cose! e~™ sin(wgt) pe 2.27 
n(t)| |e “Ly cos(wgt) + wg sin(wgt)} e "Lwgcos(wgt) — y sin(wgt)] |L C2] ed) 


Substituting Equation (2.26) into Equation (2.27), this solution can be written in terms 


of the initial conditions: 





a ae oe 
ni] 1° "L costo + gz sino) | Lee sino 10) 5 6 
nt) | en “ Pe n(o) | O7” 
——2¢ ” sin(wt) e cos(w,t) Dy sin(@ 4f) 
d 
Letting 
nt) 
Z, — . 2.29 
A?) bad ie 
and 


= Le 
e "| cos(w gl) + on sin(w) | a eo sin(@ Jf) 





6a) 
oS va - a , (23g) 
= a oe sin(w gt) é | cos(w ft) Ta Oq sin(w) | 
the solution can be written: 
Zt) = O{DZ{1) (2.31) 


where Q, is the state transition matmx of the i* mode. The nonhomogeneous solution 
is obtained from: 


Z(t) = O(NZA2) + Vx; FO (2.32) 


where the discrete-time input matrix 1s given by: 
7 
r= | ®(2)Bér, (2.33) 
0 


B = (0 1]" is the input matrix for the continuous-time system, and T is the sampling 
time [Ref. 4: p. 59]. Solving Equation (2.33) vields: 





L E —€))' costa i = eT sin(wg) | 
P= leee , (2.34) 


| -yT.. 
ra " sin(wyl) 


The discrete-time state equation for the 7* equation of motion can now be written: 
Z(kT + 1) = O(T) Z{kT) + 1 (T)x; FT) (2.35) 


where @, and I’, are evaluated at t= T. Summarizing the terms: 
e Z,1s a vector of the i* modal amplitude and the /* modal velocity 
e @, is the 7 state transition matrix 
e | ,is the é* input vector 
e xi is the transpose of the i* mode shape vector 
e F 1s the contro! torque force vector applied at a point 


e 7 1s the sampling time 


e 6k ys the time index. 


Equation (2.35) can be expanded to include a disturbance input: 
ZA{KT + 1) = O(T) Z(KT) + F(T)x; CF(KT) + WkT)] (2.36) 


where w(kT) 1s the disturbance input. Equation (2.36) is the discrete-time mathematical 
modal model describing the motion of the space structure in terms of its natural modes 
of vibration. The computer simulation of Chapter IV will solve the equations of motion 


by iterating this discrete-time model. 


HI. THE CONTROL SOLUTION 


A. INTRODUCTION 

The mathematical model derived in Chapter II represents a system with a control 
input and a random disturbance input. The problem that must be solved now is the 
development of the input that will effectively control the structure. The plant order is 
sufficiently large to preclude all but an optimal approach to the control problem 
[Ref. 4: p. 337]. Optimal control techniques attempt to find a control law which forces 
the system to follow a path that minimizes a given performance measure. [Ref. 5: p. 11] 

Control system optimization involves the selection of a performance measure that 
describes a given characteristic, or property, of the system. Minimization of the per- 
formance measure is achieved by calculating the control gains and applying them to the 
system through state feedback. A suitable performance function based on the 
vibrational energy of the structure will be developed in Section B. This performance 
function, also referred to as the cost function, will provide a means of comparing various 
control models. The solution to the discrete-time Riccati equation, providing the nec- 
essary feedback gains required for an optimal control scheme, will be presented in Sec- 
tion C. A control system for the space station is then generated bv applving the optimal 


control solution developed in this chapter to either the full order or reduced order model. 


B. SYSTEM PERFORMANCE EVALUATION 
The system mathematical model has been established. An expression must now be 

derived that will quantify control system performance and allow for a comparisons be- 
tween competing systems. Kirk states: 

In selecting a performance measure the designer attempts to define a mathematical 

expression which when minimized indicates that the system is performing in the 

most desirable manner. Thus, choosing a performance measure js a translation of 

the systems physical requirements into mathematical terms [Ref. 5: p. 34]. 
Structural vibration is a physical property of the space station that disrupts the envi- 
ronment within the system. Vibrational energy in the structure can be due to any 
number of disturbances, for example: 

@ space shuttle docking 

e rotating machinery 


e positioning jet reaction 
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® movement of crew members 


® rotation of solar panels. 


A reasonable performance function is based on the vibrational energy of the structure. 
This performance function results in the control forces needed to minimize this vi- 
bration. 


The performance measure 1s placed in a precise mathematical framework: 
J = E{T.E.] + E[u’Ru] (3.1) 


where 
e E[ + ] is the expected value of « 
e T.E. is the total vibrational energy 


e u'Ru limits the magnitude of the control force. 


The total energy consists of the potential and kinetic energies of the structure at any 


point in time. This can be expressed in matrix form: 
= =e. T AE 
T.E() = P.E() + K.E() = > [qa ()Kq() + 4 (DMQ(9 J. (3.2) 


Applying the coordinate transformation of Chapter II, the energy equations are written 


in terms of the modal amplitudes and modal velocities as: 


P.E(K) = +) waits) (3.3) 
i=] 
and 
K.E.(k) = > 72(k). (3.4) 
i=] 


Note that the energy terms are summed over the m modes of the system. Defining the 


State weighting matrix Q, as: 


2 
Q = Ke 7 (3.5) 


permits the total energy to be written in terms of the modal state vectors as: 


Ll 


: i(k) 
T.E(K) = >) [nl nk] Q ” oy (3.6) 
i=] i 


and applying Equation (2.29) results in: 


T.E{k) = +) ZIWQ Z,(k). (3.7) 
i=] 


The cost function of Equation (3.1) can now be written: 


J= E| = » Z; (k)QZ{k) + uF ORuCA | (3.8) 
i=] 


where J is a quadratic performance measure. A great deal of theory is available con- 
cerning the solution of quadratic optimal control problems. [Ref. 6: p. 84] 

The system performance will be found by computer simulation and comparisons of 
various reduced order models will be made. Evaluation of the expected values by monte 
carlo methods is tedious. An alternative approach is to compute the expected values 


from the impulse response: 


(oe) 1 4! 
J = >| > aCe + hT (kh, (k)] + HECORHAC | (3.9) 


where 
e h,(k) 1s the impulse response of 7 
¢ h,(k) is the impulse response of 1 


e h,(k) is the impulse response of u. 
This equation can be easily evaluated by computer simulation. [{Ref. 6: p. 85] 


C. RICCATI SOLUTION 

The solution to the quadratic optimal control problem 1s well known and extensively 
documented [Ref. 4: p. 338]. The optimal gain matrix must be chosen so Equation (3.8) 
will be minimized when the feedback loop is closed in Equation (2.36). The optimal 
control is state feedback: 


F(k) = u(k) = LZ(k) (3.10) 


where the optimal gain matrix, L, 1s computed: 
are eee) SLX ne (EX) S@ (3.11) 


and 
e S is the steady-state matrix solution of the Riccati equation 


e XI is the mode shape matrix of the control node. 


The product of the gain matrix and the time varying state matrix, Z, defines the contro! 
torque vector, u(k). 
The solution matrix, S, is found by solving the discrete-time Riccati equation!: 


S=0'S® — O'SPX)(R + (TX))'SPX!)'(rx!)'so + Q. (3.12) 


The solution to the Riccati equation is achieved by a variation of the Hamultonian- 
Eigenvector approach. [Ref. 7,8] 

The physical control of the structure 1s achieved by a system of control moment 
gyros or torquers. The mode shape vector X,; comprises the six degrees of freedom of 
the :* mode: three degrees of modal deflection and three degrees of modal slope. For 
these point torquers, only the modal slopes at the torquer location are of concern in the 


development of the control input vector. The input matrix can then be simplified: 
B= TX) =x’, (3.13) 


Where x’, is the row matrix of the modal slopes at the control moment gyro location. 


This simplification can also be applied to the point input disturbance where: 
BN = IX’, (3.14) 


and 
e BN is the noise input matrix 


e x’, 1s the modal slope matrix associated with the disturbance node. 


The discrete-time state of Equation (2.36) can now be written as: 


1 The subroutine to solve the discrete-time Riccati equation was provided by Prof. Alan J. 
Laub, University of California, Santa Barbara. 


Z(k +1) = OZ(k) + Bu(k) + BN w(k) (3.15) 


where w(k) is the applied disturbance input. 

The space station has now been completely defined. It exists as a discrete-time 
mathematical model, consisting of a set of /¢ uncoupled simultaneous differential 
equations of motion, expressed in terms of the natural modes of vibration. The control 
svstem is defined by the performance measure and either the full order or reduced order 
model. The vibrational energy of the structure is the measure of performance by which 
the reduced order control system will be judged. Input to the discrete-time Riccati 
equation consists of the state transition and input matrices of the model as well as the 
state weighting matnx and the modal slopes. The optimal control torques are based on 
the solution of this equation. The final product is Equation (3.15) which forms the 
framework for the simulation of Chapter IV. 


IV. SIMULATION AND RESULTS 


A. INTRODUCTION 

The objective of the simulation 1s to determine the system response to disturbances 
applied at various points on the structure. This is accomplished by iterating the 
discrete-time model developed in Chapter II. 

The process outlined above is covered in three sections. The specifics of the data 
used to model the space structure are discussed in Section B. An overview of the simu- 
lation program is given in Section C, and the specific disturbance used to excite the 
structure is introduced. The results of the simulation for various control conditions and 


disturbance locations are presented and discussed in Section D. 


B. MODEL DATA 

The dynamic model used in the simulation is for a preliminary space station config- 
uration; the phase II dual keel structure2. The model consists of the first 100 natural 
modes, i.e., the natural frequencies for the first 100 modes and a matrix of 100 mode 
shapes for 114 nodes with six degrees of freedom at each node. The natural frequencies 
“of the structure are shown in Table 1. The first six natural frequencies correspond to 
rigid bodv modes with the first bending mode beginning at number seven. Development 


of the reduced order control will be based on these bending modes. The modal ampli- 
lb-s? 
in 





tude data was normalized for a modal mass of 1 , specifying that the units for the 
simulation be given in the English svstem. 

The svstem of 114 nodes is quite large, and 1t would be difficult to observe the re- 
sponse due to disturbances at each node. Therefore, three nodes of particular interest 
are considered in the simulation. They are the shuttle docking point at node 23, the 
alpha-joint3 at node 55, and the location of the control moment gyros at node 69. The 
relative locations of these nodes on the structure is depicted in Figure 2 on page 17. 
These particular nodes are singled out for consideration because they are the point lo- 


cations of either the applied control forces (node 69) or the applied disturbance force 


2 The dynamic model for a preliminary space station configuration is provided courtesy of 
McDonnell Douglas Astronautics Company, 5301 Bolsa Avenue, Huntington Beach, CA 92647. 


3 The alpha-joint is the physical connection between the fixed structure and the rotating solar 
panels. 


(nodes 23 or 55). The modal slopes associated with these particular nodes are listed in 


Appendix B. 


Table 1. 
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C. SIMULATION PROGRAM 

The Space Structure Simulation Program, listed in Appendix A, is written to simu- 
late the space structure using the discrete mathematical model developed in Chapter I] 
and the applied control developed in Chapter III. It provides the iterative solutions to 
Equation (3.15) and computes the cost function, Equation (3.8), using the equivalent 
form, Equation (3.9). 

The program is structured in block format with each block providing the major in- 
put to the next. The major blocks, or sections, are: 

¢ screen interaction 


e establishment of the required matrices 


solution of the discrete-tume Riccati equation and calculation of the optimal control 
gain matrix 


time iteration. 


A discussion of these blocks will follow. 





Figure 2. Node Location Diagram. 


The simulation program begins with an interactive phase. This capability allows the 


user to select the values of certain parameters that are not fixed by either the program 


or the input data. The parameters that may be varied are: 


maximum size of the modal model 

number of control modes to be used, 1.e., the size of the reduced order model 
node location of the applied disturbance 

axial direction of the applied disturbance 


sampling time 
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¢ damping factor 


e value of the diagonal elements, r, of the control weighting matrix, R. 


The choice of values for these items will fully specify the entire simulation. Several of 


these values were kept constant in generating the results presented in this thesis. They 


are: 
® initial value ofr: ehGse 
¢ sampling time: 1x10-2 
e damping factor: 1x10-3. 


The value of r was picked based on trial and error, and represents a compromise between 
the control cost and the vibrational energy of the structure. The sampling time was 
based on the period of the highest frequency of the system. The chosen value represents 
a sampling frequency that is approximately ten times faster than the highest natural 
frequency, resulting in a minimal amount of aliasing. The damping factor was chosen 
to yield a lightly damped structure. This choice is based on previous space station 
studies that have used values in the range of 0.0001 to 0.005 [Ref. 9,10]. 

The program proceeds to construct the matrices of Equation (3.15) based on the 
choice of parameters, the natural frequencies of Table 1 on page 16, and the modal 
slopes at the control node and the disturbance node. The state weighting matrix, Q, is 
formed in accordance with Equation (3.5), and the control weighting matrix, R, 1s de- 


termined by the users selection of an appropriate value r where: 


100 
R=r/0 1 Ol. (4.1) 
001 


The next section of the program has two functions. The first is to obtain the 
discrete-time Riccati solution [Ref. 7,8] as given by Equation (3.12). This is accom- 
plished by applying the subroutine, RICDSD, to the appropriate matrices of the previ- 
ous paragraph. The second function is to apply the solution to obtain the optimal 
control gain matrix. A sample of this portion of the program is shown in Figure 3 on 
page 19. This sample includes the Riccati solution matrix, the closed loop eigenvalues 
of the reduced order model, and the control gain matrix. This information is provided 
as an output file for any sized reduced order model. At this point, all constant value 


matrices have been established. 


STARTING MODE NUMBER: 

NUMBER OF MODES SCANNED: 

LAST CONTROLLED MODE: 

NOISE INPUT NODE: 5 

INITIAL R VALUE: 0. 1000E-11 
SAMPLING TIME: 0. LOOOE-01 
DAMPING FACTOR: O21000E =02 
OBSERVATION TIME: 120.0 MINUTES 
SIZE OF MODAL MODEL: 100 MODES 


THE RICCATI SOLUTION IS: 
0.593595E+02 0.402677E+01 ~. 207784E-01 1395598, =O 1 
0.402677E+01 0. 817013E+01 - 359488E-01 «. 728645E-01 


-. 207784E-01 -. 359488E-01 0. 619850E+02 0. 485536E+01 
-. 349359E-01 =9/25045E5-01 0. 485536E+01 OS 95792Z9E+01 


THE CLOSED LOOP EIGENVALUES ARE: 


Real Imaginary 


O67 555600175230962835 . DCOOODD000000000000E+00 
O769377679206653312 . DCOODODD000000000000E+00 
0. 994092984040947231 . HOODDDD00000000000E+00 
0. 994306266098994654 . DHOODD0000000000000E+00 


THE CONTROL GAIN MATRIX L IS: 
se OZo1SETUG wevo29 5 1E+0G -. 647580E+05 pee d27 boo 


Oes2127 95705 0.617015E+05 -. 109143E+06 -. 205384E+06 
O77 1026E+05 O. 110078E+06 =. 509676E+06 22959551 EUG 





Figure 3. Sample Output for a Two Mode Model of Reduced Order Control. 


The time iterating portion of the simulation comprises the final major block of the 
simulation. It consists of a number of sub-sections designed to accomplish a specific 
sequence of tasks. The first sub-section computes the control torques based on the 
control or feedback gain matrix and the time varying state vector. The control torques 
form a 3x1 vector and are obtained by the process of Equation (3.10). For the initial 


conditions at k= 0, the control is: 


0 
u(0) = LZ(0) = | 0}. (4.2) 
0 
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A new torque vector will be computed for each increment of k. 

The disturbance input used in the simulation 1s a unit impulse. The impulse imparts 
a disturbance, at time k=0, to a selected point (node 23, 55, or 69) on the structure, and 
the system response is then observed for the desired period of time. The disturbance 
may be applied in the x, y, or z direction, depending on the users choice during screen 
interaction. Provisions are also made within the program to obtain an uncontrolled 
impulse response. 

The system cost function comprises the next sub-section of this block. The cost 
function is computed using Equation (3.9) with the total energy (T.E.) of the structure 
summed over all i modes for each k* increment of time. This energy is combined with 
the control cost, u’Ru, to vield the cost function, J. In addition to the system cost, this 
section computes the total energy for each mode. The graphical analysis in Section D 
is based on the tabulated results of this sub-section. 

The final sub-section of the block is the update of the discrete state equation. 
Equation (3.15) 1s updated over i modes for each k” time increment and the new state 


vector Z results. For time k=O, this equation can be written: 
Z(1) = ®Z(0) + Bu(O) + BN 6(0) (4.3) 
where 6(0) 1s the applied unit impulse. When k > 0, Equation (4.3) reduces to: 
Z({k +1) = OM Z(k) + B u(k). (4.4) 


The program loops back to the beginning of the time-varying block when the update of 
the state vector is complete. A new torque vector 1s computed, and the process de- 
scribed above continues until the response of the system has been found over the desired 
ume period. 

When the simulation is complete, the user is prompted to choose between running 
again with a new control weighting matrix, R, or termination of the simulation. The 
choice of a new weighting matrix reinitiates the simulation process with all parameters 
remaining the same except for the matrix R. 

Output data can be tailored to fit user requirements. The program 1s written to 
provide the total system energy, the control cost, and the total energy per mode for a 


control based on a given reduced order model. 
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D. RESULTS | 

The response of the svstem to a disturbing force is provided by the simulation of 
Section C. The simulation is used to compute the systems response to a disturbance 
applied at either the space shuttle docking point, node 23, or the alpha-joint, node 55. 
These points frequently experience large disturbance inputs. 

The simulation provides the energy of the vibrating structure in two formats. First, 
the total energy of the structure, as defined by Equation (3.8), is obtained. This provides 
a single value for a given condition of control. Second, the total energy of each mode 
is Obtained. Essentially, this breaks down the single total energy value into 100 com- 
ponents, each representing the total energy of an individual mode. Tabulating the en- 
ergy per mode will permit a more detailed look at the system response and the transfer 
of energy between the control system and the structure. Figure 4 on page 22 shows the 
full order system response due to a disturbance applied at each of the nodes of interest. 
Observe that the cost of control in terms of energy increases as the size of the reduced 
order control increases, but then begins to decrease. The important point to note is the 
size of the reduced order control required to bring the cost back to, and below, the un- 
controlled level. 

A major drawback arose during the simulation process. To obtain a single data 
point of Figure 4 on page 22 it required 40 minutes of computer CPU time to simulate 
the system. If the size of the model could be reduced and a similar response shape ob- 
tained, then the svstem could be simulated more efficiently. 

Because high frequencies damp out quicker, the system was truncated at 56 modes; 
reasoning that the higher modes will die out quickly and contribute little to the cost. 
This size model would allow for up to 50 modes of control. Combined with the s1x rigid 
body modes, the system could then be simulated for full control, something that could 
not be achieved with a full order model. The model was truncated, and the results are 
shown in Figure 5 on page 23. The response very closely matches the shape of 
Figure 4 on page 22 for less than 30 modes of control. Above 30 modes, the comparison 
is less exact, but 1s still decreasing in energy. The significance of the truncated model is 
that it requires one-fifth of the CPU time (eight minutes vs. 40 minutes) to simulate the 
system for a two hour period. A two hour simulation period is used because this allows 
observation of the system for two time constants of the lowest natural frequency. 

To use the truncated model, the energy content of the higher modes must be negli- 


gible. A comparison of the energy plots of both models indicate that the total energy 
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Figure 5. System Energy of a Reduced Order Model with Reduced Order Control. 


of both models is nearly zero at 50 modes of control. The full order system still has 44 
modes contributing to the energy cost while the reduced order model has none. The 
conclusion is reached that a good approximation of system response can be made by 
reducing the order of the model. The savings of computer costs is considerable when a 
number of control conditions are simulated. 

A limitation of the total energy depiction of Figure 4 on page 22 and Figure 5 on 
page 23 1s that an insight to the behavior of each mode js not provided. Figure 6 on 
page 25 thru Figure 15 on page 34 provide a graphical indication of mode response to 
an impulse applied at node 55. Each figure is a depiction of the energy in each mode for 
a given size of reduced order control compared to the uncontrolled system. 

These ten figures show that the response of a directly controlled mode is essentially 
zero. They also show that as the number of controlled modes changes, there are certain 
groupings of modes that significantly increase their contribution to the cost over the 
uncontrolled system. These trouble modes appear in the areas of modes 30 to 32, 50 to 
52, and 54 to 56. All other modes show little response to a change in the number of 
controlled modes. 

A look at the system response to a disturbance applied at node 23 is shown in 
Figure 16 on page 35 thru Figure 25 on page 44. Note that the trouble modes of a node 
23 response are the same as for a node 55 response. Although the shape of the uncon- 
trolled and the controlled responses differ between the nodes, the offending modes re- 
main the same. 

The results indicate that there are certain modes, or groups of modes, that are 
troublesome. The same modes cause problems irregardless of whether the disturbance 
is applied to node 23 or node 55. The trouble modes must therefore be modes with large 
coupling to the control torquers. Two things influence this coupling, the natural fre- 
quency of the mode and the modal slope at the torquer location. As the natural fre- 
quency goes up the damping goes up proportionally but the spread between frequencies 
is small so damping is comparable for all modes. The modal slopes for the trouble 
modes (see Appendix B) indicate larger magnitude than those of the surrounding modes. 
The effect of these slopes is to increase the influence factor of the torquers. This can 
lead to a large excitation of the modes by the control system. If this is the case, it may 
be possible to identify possible trouble modes by comparing natural frequency spread 
with control node slope magnitude. A combination of a small frequency spread and a 


large slope magnitude could tag potential problem modes before simulation. 
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Figure 6. Node 55 Response with 5 Controlled Modes. 
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Figure 7. Node 55 Response with 10 Controlled Modes. 


26 


Node 55 Response 
Energy Per Mode 


ome 
a 
LEGEND : 
° Uncontrolled 

AS controlled modes 


e 
8 
% 
% 
8 
es 


Energy (in—lbs) 
15.0 200 250 300 


10.0 


9.0 


0.0 





—5.0 


0.0 10.0 20.0 30.0 40.0 50.0 60.0 
Mode 


Figure 8. Node 55 Response with 15 Controlled Modes. 
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Figure 9. Node 55 Response with 20 Controlled Modes. 
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Figure 10. Node 55 Response with 25 Controlled Modes. 
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Figure 11. Node §5 Response with 30 Controlled Modes. 
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Figure 12. Node $5 Response with 35 Controlled Modes. 
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Figure 13. Node 55 Response with 40 Controlled Modes. 
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Figure 14. Node 55 Response with 45 Controlled Modes. 
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Figure 15. Node 55 Response with 50 Controlled Modes. 
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Figure 16. Node 23 Response with § Controlled Modes. 


35 


Node <3 Response 
Energy Per Mode 


i ee 
7 @ 
LEGEND 
© Uncontrolled 
é controlled modes 


5.0 6.0 


Energy (in—lbs) 
2.0 3.0 4.0 


1.0 


0.0 





=1 0 


0.0 10.0 20.0 30.0 40.0 50.0 60.0 
Mode 


¢ 


Figure 17. Node 23 Response with 10 Controlled Modes. 
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Figure 18. Node 23 Response with 15 Controlled Modes. 
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Figure 19. Node 23 Response with 20 Controlled Modes. 
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Figure 20. Node 23 Response with 25 Controlled Modes. 
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Figure 21. Node 23 Response with 30 Controlled Modes. 
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Figure 22. Node 23 Response with 35 Controlled Modes. 
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Figure 23. Node 23 Response with 40 Controlled Modes. 
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Figure 24. Node 23 Response with 45 Controlled Modes. 


43 


Node <3 Response 
Energy Per Mode 


*10° 
8.0 


LEGEND 
Uncontrolled 


7.0 


controllea Mmoaes 


6.0 


5.0 


4.0 


Energy (in—lbs) 
3.0 


2.0 


1.0 


0.0 


—1.0 





0.0 10.0 20.0 30.0 40.0 50.0 60.0 
Mode’ 


Figure 25. Node 23 Response with 50 Controlled Modes. 


The results, total system energy and total energy per mode, have been presented in 
a graphical format. Total system energy provides a good indication of total perform- 
ance, but not much more. Figure 4 on page 22 indicates a large number of modes must 
be controlled to reduce the energy cost to a level below that of an uncontrolled system. 
An alternative approach to the control problem would include the control of the trouble 
modes in the formulation of the cost function, which will result in control gains which 
do not excite these modes. 

A better means of observing the system is achieved by presenting the results in the 
alternate format of energy per mode. Trouble modes identified themselves in Figure 6 
on page 25 thru Figure 25 on page 44. These groups of modes were sensitive to any 
change in the control structure of the space station, and are major contnibutors to the 
energy cost until they become directly controlled. Reaching a desirable energy cost will 
require a reduced order control system of at least 50 modes. This may be prohibitive 


with regard to complexity, cost and weight. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

A mathematical model of a large space station was developed, several optimal con- 
trol solutions based on reduced order models were found, and the effects of reduced or- 
der modeling on the control of a large space structure have been presented in this thesis. 
It was shown that a straightforward approach to reduced order control can be achieved 
for a space station. The negative aspect of this approach is that the number of modes 
required to reach an acceptable level of control is considered to be excessive from a cost 
and weight viewpoint. 

The existence of certain trouble modes was observed. These mode groups contrib- 
uted significantly to the cost of control until the control model became large enough to 
directly control these modes. Also, 1t was demonstrated that these mode groups re- 


mained troublesome regardless of the location of the disturbance. 


B. RECOMMENDATIONS 
Based on the results of this work, there are several additional areas that mav be in- 
vestigated with respect to the effects of reduced order control on a space structure. 


Some of these areas are: 


e to develop a cost function that considers control of the trouble modes in the opti- 
mal gain solution 


® to examine the effects of notch filtering to supress excitation of the trouble modes 


® to utlize Karhunen-Loeve expansion methods for reduced order control modeling. 


Research in these areas may produce results more favorable to solving the control 


problems of a large space structure. 
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APPENDIX A. SPACE STRUCTURE SIMULATION PROGRAM 


LOAVAACA PAT ACTA ACACH TCR ICA ACHE PACHTER AICTE AHA AKA KAR aRICK KR RRR 


Wietick eke 
AT Bas SPACE STRUCTURE SIMULATION PROGRAM Feieieieye 
tecicit tetctick 


POPC ACA TACHA FCA AACE Re FCPCACR ERR RRR AAAI RARER 


FATT AACA AACA ACHE ACA CHERCTCR TCA CH REE ACTCT TAC PPA ICH ACA ATACHK ACHAEA RICK 


Feiciriek VARIABLE DECLARATIONS EKKK 
PRAHA AIAN IAI AAAI AAA IA IARI IAA AAAI AA AAA AAA AV AAAAARAK 


EXTERNAL EXCMS,RICDSD 

CHARACTER*6 NAM 

CHARACTER*1 AGAIN, CORECT,RAGAIN 

INTEGER NODE,MODE,KQ,EMODE , SMODE ,CTADJ,CT,CF,KA, I 
INTEGER LOOP, PRNT,MODAL, COUNT, PRNTG,NF,NG,NH,NZ,V 
INTEGER IPVS(100) , ITYPE( 200) 


heen 


REAL RTOTAL,RMODEN(7: 100) 

REAL*8 PHII(2,2,100),GAMMA(2,100),EGT,GMA,WN,W1,X1T,X2T 
REAL*8 PHI( 188,188) ,B(188,3),BN(188,3),R(3, 3) ,RR(3, 3) 
REAL*8 RINV(3,3),RRINV(3,3),X1(7: 100) ,X2(7: 100) ,MODEN( 7: 100) 
REAL*8 COSW1T,SINW1T,COST,CNTCST, ENERGY , TOTCST, RM 
REAL*8 TCX,TCY,TCZ,DAMP,SAMPT, PI,SUM1,SUM2,SUM3,SUMC 
REAL*8 TNX,TNY,TNZ, IMPX,IMPY, IMPZ,IMPLSX, IMPLSY, IMPLSZ 
REAL LAMA(100) , UGVEX(100,3),RNODE,RMODE,MIN, TIME, SAMPTM 
REAL UG69(100,3) ,UG23(100,3) ,UG55(100,3) 

REAL*8 H(100,100),G(100,100) ,L(3,100) ,BT(3, 100) 

REAL*8 Z(200,200) ,W(200, 200) ,ER( 200) ,F(100, 100) ,EI(200) 
REAL*8 SCALE(200),TEMP(100,3),TEMP1(3, 100) ,WORK( 100) 


Prete RHEE ERE RR RE RE RRR RER REE RIR EIR EERIE TRE 


Wickit VARIABLE DEFINITIONS Wrecker 
derederdetodedeve de dete Fee RICK IK TKI IEE IK UIA NIK III IPI IAAI REI 


LAMA = VECTOR OF THE SQUARE OF THE NATURAL FREQUENCIES 

UGVEX = MODE POSITONS AND SLOPES OF THE NODAL POINTS 

PHII = TRANSITION MATRICIES FOR EACH MODE 

PHI = BLOCK DIAGONAL STATE TRANSITION MATRIX CONSISTING OF 
THE INDIVIDUAL PHII MATRICIES 

GAMMA = INPUT TRANSITION MATRIX 

B = INPUT MATRIX OF GAMMA AND CONTROL NODE SLOPES 

BN = NOISE INPUT MATRIX OF GAMMA AND NOISE NODE SLOPES 

DAMP = DAMPING FACTOR 

SAMPT = SAMPLING TIME 

IMPLSE = IMPULSE INPUT FUNCTION 

TCX, TCY, TCZ = CONTROL TORQUE VALUES 

IMPX, IMPY, IMPZ = AXIS IMPULSE NOISE VALUES 

ENERGY = SYSTEM ENERGY COST VALUE FOR A GIVEN POINT IN TIME 

CNICST = SYSTEM CONTROL COST VALUE FOR A GIVEN POINT IN TIME 
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Cie 


oe he) 


7) C2 c2 © 


COST = TOTAL SYSTEM COST VALUE FOR A GIVEN POINT IN TIME 
TOTCST = SYSTEM COST SUMMED OVER ALL TIME 
MIN = NUMBER OF MINUTES SYSTEM WILL BE OBSERVED 


RIKI HK SAMPLE OF SPACE EXEC FILE KicteK Ici KK 


THIS FILE MUST BEGIN IN COLUMN 1 AND RUN IN THE FOLLOWING 
SEQUENCE FOR THE INITIAL RUN OF THE PROGRAM: 


FORTVS SPACE (COMPILES PROGRAM) 
SPACE (EXECUTES EXEC FILE) 
LOAD SPACE (START (LOADS AND EXECUTES PROGRAM) 


SUBSEQUENT PROGRAM RUNS CAN ELIMINATE "FORTVS SPACE" IF NO 
CHANGES HAVE BEEN MADE TO THE PROGRAM AND CAN ELIMINATE 
RUNNING THE EXEC FILE. 


FI 4 DISK THESIS INPUT A (PERM 

FI 30 DISK X1 OUTPUT A (RECFM F BLOCK 80 PERM 

FI 31 DISK MOLENG OUTPUT A (RECFM F BLOCK 80 PERM 
FI 32 DISK TORQUE OUTPUT A (RECFM F BLOCK 80 PERM 
FI 33 DISK ENERGY OUTPUT A (RECFM F BLOCK 80 PERM 
FI 34 DISK MDECST OUTPUT A (RECFM F BLOCK 80 PERM 
FI 35 DISK COUNT OUTPUT A (RECFM F BLOCK 80 PERM 
FI 40 DISK UTILITY OUTPUT A (RECFM F BLOCK 80 PERM 
FI 41 DISK RUN OUTPUT A (RECFM F BLOCK 80 PERM 

FI 42 DISK NODE69 INPUT A (RECFM F BLOCK 80 PERM 
FI 43 DISK NODE23 INPUT A (RECFM F BLOCK 80 PERM 
FI 44 DISK NODE55 INPUT A (RECFM F BLOCK 80 PERM 


THE THESIS INPUT FILE CONTAIN THE NATURAL FREQUENCIES OF THE 
SYSTEM (LAMA VECTOR). MODAL SLOPES FOR THE DESIRED NODES ARE 
CONTAINED IN THE INPUT FILES NODE69, NODE55 AND NODE23. THE 
REMAINING FILES IN THE EXEC CAN BE STRUCTURED TO HANDLE OUTPUT 
AS DESIRED. 


FRR EERE KEKE ERK ERE RERKREREREREEEREERERRRIEKE 


PI = 4.0D0 * ATAN(1. ODO) 


SAMPT = 0.0 
DAMP = 0.0 
MODAL = 0 
NF = 100 
NG = 100 
NH = 100 
NZ = 200 


week =NUMBER OF MINUTES THE SYSTEM WILL BE OBSERVED **%#**r* 
MIN = 120.0 


FIAT AAAI IATA IAAI IASI AA AAAAAR AAAS AADAC AISA AAI AAAI 
ToanaiaNre SET LENGTH OF MODAL MODEL ae 
FAA A AHIR AAA IAT IIT IAAI AAAI IAA A AIA AAAI AAAI 


CALL EXCMS ('CLRSCRN' ) 
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WRITE (6,1008) 
WRITE (6,*) ' SET MAXIMUM OPERATING MODE NUMBER’ 
READ *, MODAL 


ce 
C Sesetevedede teve tetete tevtete teste fede Fe T TT TCII TRTIC T TT ICH RRRERRRRERRIRR RRR 
C Seve devede READ LAMA MATRIX teteve steve 
O) Sea scdetede detested desetedeste kote Keele Fesete sede te FINI IKE IKKE KI REIKI 
C 

CALL EXCMS ('CLRSCRN' ) 
c 

READ(4,1001) NAM 

READ(4,1002)( LAMA(I) , I=1, 100) 
6 
C POPPA ACH AACA CE TP ORE HEPAT ACA ACHAEA TAR IARRAI ERR EKEKERE RC RRR RRR aE 
G deiekiek SCREEN INTERACTION we 
c FEV HE FOP HAHA HCH CH FETA TCA PTCA CSE FPA ICTRACKER RK RR KKK 
C 
500 CALL EXCMS ('CLRSCRN' ) 
C 
‘& tovededtetedetetedeve seve STARTING MODE NUMBER tekecdkkKkickiKk 
c 


WRITE (6,1004) 


10 WRITE (6,*) 'ENTER THE STARTING MODE NUMBER(MUST BE 7 OR GREATER 
oe 
READ *, RMODE 
IF(MOD( RMODE,1.0).NE.0) THEN 
WRITE (6,*)' MODE SELECTION MUST BE AN INTEGER VALUE. RE, 
+ENTER’ 
GOTO 10 
ENDIF 
SMODE = INT( RMODE) 
IF (CSMODE. LT. 7). OR. (SMODE. GT. 100) )THEN 
WRITE (6,*)'MODE CHOICES ARE LIMITED TO 7 THRU 100!! REENTER. ' 
GOTO 10 
ENDIF 
C 
C HAEKKKKKK TE NUMBER OF MODES USED FOR CONTROL HIRE ET 
C 
15 WRITE (6,*) ‘ENTER THE NUMBER OF MODES FOR CONTROL. THE CHOICE MU 


+ST BE FROM 1 TO 94: 
READ *, RMODE 
IF(MOD(RMODE ,1.0).NE.0) THEN 
WRITE (6,*)' MODE SELECTION MUST BE AN INTEGER VALUE. RE 
+NTER'’ 
GOTO 15 
ENDIF 
MODE = INTC RMODE) 
IF (CMODE. LT. 1).OR. (MODE. GT. 94) )THEN 
WRITE (6,*)'MODE CHOICES ARE LIMITED TO 1 THRU 94!! REENTER. ' 
GOTO 15 
ENDIF 
EMODE = SMODE + MODE - 1 


RR RRR NOISE INPUT POSITION PARTITE IRAE 


qaqa 


CALL EXCMS ('CLRSCRN' ) 
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qaqa 


maa ca 


wooo 


WRITE (6,1004) 


IMPX = OQ. ODO 
IMPY = 0. 0DO 
IMPZ = 0. 0DO 


WRITE (6,*) ‘ENTER THE NOISE INPUT LOCATION. 


+55 OR 69. 
READ *, RNODE 


IF(MOD(RNODE,1.0).NE.0) THEN 
WRITE (6,*)' NODE SELECTION MUST BE AN INTEGER VALUE. RE 


+NTER' 
GOTO 20 
ENDIF 
NODE = INTC RNODE) 


IF (( NODE. LT. 1). OR. (NODE. GT. 114) ) THEN 


THE CHOICES ARE 23, 


WRITE (6,*)'NODE CHOICES ARE LIMITED TO 1 THRU 114!! REENTER. ' 
GOTO 20 

ENDIF 

WRITE (6,*) 'SELECT THE NUMBER OF THE NOISE INPUT AXIS ' 

WRITE (6,*) ‘ ' 

WRITE (6,*) ' 0 NO NOISE INPUT' 

WRITE (6,*) ' 1 X AXIS INPUT ' 

WRITE (6,*) ' 2 Y AXIS INPUT ' 

WRITE (6,*) ' 3 Z AXIS INPUT ' 

WRITE (6,*) ' 4 INPUT ON ALL AXIES ' 

READ *, AXIS 

RRR RRR RI R MATRIX VALUE Tietietetesen RRR eR 

CALL EXCMS ('CLRSCRN' ) 

WRITE (6,1004) 

WRITE (6,*) 'ENTER YOUR INITIAL R VALUE: ' 

READ *, RM 

Sotoseseyvesesesesensese se SAMPLING TIME Sededertetedesetededeiciese 


WRITE (6 ,*) ‘ENTER SAMPLING TIME (MUST BE EQUAL TO OR LESS THAN 0. 


+04 SEC): 
READ *, SAMPT 


SAMPTM = ((2. ODO*PI) /SQRT( LAMA(100)))/2. ODO 
IF (SAMPT. GE. SAMPIM) THEN 
ERROR: SAMPLING TIME MUST BE LESS THAN OR E 


WRITE(6,*)' 
+QUAL TO' 
WRITE(6,1005) 
GOTO 25 
ENDIF 


Rik KK KT 


SAMPTM 


DAMPING FACTOR 


WRITE (6,*) ‘ENTER DAMPING FACTOR D: ' 


READ *, DAMP 


IF ( (DAMP. LT. 0). OR. (DAMP. GT. 1. 0) ) THEN 
WRITE (6,*)'DAMPING FACTOR IS LIMITED TO 0 THRU 1.0 - 


GOTO 30 
ENDIF 
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KEKEKKEKEEREKREE 


REENTER' 


520 


706 
707 
708 
709 
710 
“let 


KiIewseteseddede tctr VALUES CORRECT 7??? MRietetesdetededededesedeser se 
CALL EXCMS ('CLRSCRN’ ) 

WRITE (6,1004) 

WRITE (6,*) ‘ARE THESE VALUES CORRECT? (Y/N) 

WRITE (6,*) ' ' 

WRITE (6,706) SMODE 

WRITE (6,707) MODE 

WRITE (6,708) NODE 

WRITE (6,709) RM 

WRITE (6,710) SAMPT 

WRITE (6,711) DAMP 

FORMAT (' ',' STARTING MODE NUMBER: ',I2) 
BORMAT.(.* , NUMBER OF MODES SCANNED: ‘',I2) 
FORMAT (' ‘',' NOISE INPUT NODE: ',I3) 

FORMAT (' ‘',' INITIAL R VALUE: ‘',E12.4) 
FORMAT (' ',' SAMPLING TIME: ' E12. 4) 
FORMAT (' ',' DAMPING FACTOR: Ele) 


READ(*,1010) CORECT 
IF(CORECT. EQ. 'Y') THEN 
GOTO 510 
ELSEIF(CORECT. EQ. 'N') THEN 
GOTO 500 
ELSE 


WRITE(6,*) ‘YOU MUST CHOOSE UPPER CASE "Y" OR "N". 


+SELECT AGAIN. ' 
GOTO 520 
ENDIF 


CALL EXCMS ('CLRSCRN' ) 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 


(41,700) 
(41,701) 
(41,712) 
(ie 702) 
(41,703) 
(41,704) 
(41,705) 
(41,713) 
(6,1008) 
Ce 


SMODE 
MODE 
EMODE 
NODE 
RM 
SAMPT 
DAMP 
MIN 


Fetctededede ies 


IF( AXIS. EQ. 1) THEN 
IMPX = 1. 0DO/SAMPT 
ELSEIF(AXIS. EQ. 2) THEN 
IMPY = 1. OD0/SAMPT 
ELSEIF( AXIS. EQ. 3) THEN 
IMPZ = 1. ODO/SAMPT 
ELSEIF( AXIS. EQ. 4) THEN 


IMPX = 1. 0D0/SAMPT 

IMPY = 1. 0D0/SAMPT 

IMPZ = 1. ODO/SAMPT 
ENDIF 
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PROGRAM RUNNING' 


NOISE AXIS INPUT AND LOCATION 


°CAPS ONLY ' 


Feretetedeteresew eK 


qQaoaca 


45 
40 


65 
60 


Uy 


aqaqnaqanaqsn 


COUNT = 0 


RIKI KR I I 


DO 40 I = 1,188 


DO 45 J = 1,188 


INITIALIZE MATRICES 


PHI(I,J) = 0.0 


CONTINUE 
CONTINUE 


DO 60 I = 1,188 


DO 65 J = 1, 


B(I,J) = 0.0 


3 


BN( I,J) = 0.0 


CONTINUE 
CONTINUE 


DO 70 K = 7,1 
X1(K) = 0. 
X2(K) = 0. 


00 
0 
0 


MODEN(K) = 0.0 
RMODEN(K) = 0.0 


CONTINUE 
DO 75 I = 1,100 


Riek cK I 


READ(42,1040) (UG69(1,K) ,K=1,3) 
READ(43,1040) (UG23(1I,K) ,K=1,3) 
READ(44,1040) (UG55(I,K) ,K=1,3) 


CONTINUE 

FIA AA HAA HAI HII IKI IK ITHAKA TRIS RRR IAHR AAAI KIRIN 
KereHses BEGIN MAIN PROGRAM wickick 
1 ay ESTABLISH PHII, PHI, B AND BN MATRICES wiekiek 


PPT CRT CIT IC CTE ICR EPRI TCRICK RICKER REE FRI TER RE RIKER ERR RIKER 


DO 600 I = SMODE,MODAL 


WN = DBLE(SQRT( LAMA(I))) 


GMA 
EGT 


DAMP*WN/2. 0 
DEXP( -GMA*SAMPT ) 


W1 = DSQRT( (WN**2) -( GMA**2) ) 
COSW1T = DCOS(W1*SAMPT) 
SINW1T = DSIN(W1*SAMPT) 


IF(WN. EQ. 0) THEN 


PHII(1,1,1) 
PHITCi eT) 
Poi 1) 
iD 2a 
GAMMA( 1,1) 
GAMMA(2,I) 

ELSE 

PHII(1,1,1) = 

PHII(1,2,I) = 

PHII(2,1,I) = 


EGT*COSW1T 
SAMPT 

0 
EGT*COSW1T 


EGT*(COSW1T + (GMA*(W1**(-1)))*SINW1T) 
(W1**( -1))*EGT*SINW1T 
- CWN**2)*(W1*%( -1))*EGT*SINW1T 
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C 
C 
600 
C 
C 


C 


C 
610 
C 


PHII(2,2,1) = EGT*(COSW1T - (GMA*(W1**(-1)))*SINW1T) 


GAMMA( 1,1) = CWN*¥**(-2))*( 1. ODO-EGT*( COSW1T+( GMA/W1)*SINWI1T) ) 
GAMMA(2,1) = (W1**(-1))*EGT*SINW1T 
ENDIF 

CONTINUE 

V=1 


DO 610 K = SMODE,MODAL 


PHICV,V) = PHII(1,1,K) 
PHEGY ;Vt1) = PHII(C1,2,K) 
EneCvtl.V¥) = PHII(2,1,K) 
PRUCVEL,Vt1) = PHIIC2,2,K) 


BCV,1) = GAMMA(1,K)*DBLE(UG69(K,1)) 

B(V,2) = GAMMA(1,K)*DBLE(UG69(K,2)) 
= GAMMA(1,K)*DBLE( UG69(K,3)) 

BC V+1,1) = GAMMA(2,K)*DBLE( UG69(K,1) 


) 
B(V+1,2) = GAMMA(2,K)*DBLE(UG69(K,2)) 
B(V+1,3) = GAMMA(2,K)*DBLE(UG69(K,3)) 
V = V+2 
CONTINUE 
DO 605 I = 1,100 
UGVEX(I,1) = 0.0 
UGVEX(I,2) = 0.0 
UGVEX(I,3) = 0.0 
IF(NODE. EQ. 23) THEN 
UGVEX(1,1) = UG23(1,1) 
UGVEX(I,2) = UG23(I,2) 
UGVEX(I,3) = UG23(I,3) 
ELSEIF(NODE. EQ. 55) THEN 
UGVEX(I,1) = UG55(I,1) 
UGVEX(1,2) = UG55(I,2) 
UGVEX(1I,3) = UG55(1,3) 
ELSEIF(NODE. EQ. 69) THEN 
UGVEX(I,1) = UG69(I,1) 
UGVEX(I,2) = UG69(I,2) 
UGVEX(I,3) = UG69(I,3) 
ENDIF 
CONTINUE 


V=1 
DO 620 K = SMODE,MODAL 


BN(V,1) = GAMMA(1,K)*DBLE( UGVEX(K,1)) 
BN(V,2) = GAMMA(1,K)*DBLE( UGVEX(K,2)) 
BN(V,3) = GAMMA(1,K)*DBLECUGVEX(K, 3) ) 


BN(V+1,1) = GAMMA(2,K)*DBLE(UGVEX(K,1)) 
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NO 
© 


C20) tee OS oe >) 
© 


25 
50 


66 
61 


89 


QO >) >) qaqa 


BNC V+1,2) = GAMMA(2,K)*DBLE( UGVEX(K,2)) 
BN(V+1,3) = GAMMA(2,K)*DBLE( UGVEX(K,3)) 
V = V+2 

CONTINUE 

CONTINUE 


KKKKKKKAKKK ESTABLISH H, F AND R MATRICES 


DO 50 I = 1,NH 
DO55J=1 
H(I,J) 
F(I,J) 
eci) 
IF(I. LE. 3) THE 
L(G) =0.6 
ENDIF 
CONTINUE 
CONTINUE 


me 


tt tt te 


DO 61 I = 1,3 
DO 66 J = 1,3 
R(1,J) = 0.0 
RINV(I,J) = 0.0 
CONTINUE 
CONTINUE 


KQ = 1 

DO 80 K = SMODE,EMODE 
HCKQ,KQ) = DBLE(LAMA(K) ) 
HC KQ+1,KQ+1) = 1. 0D0 
KQ = KQt2 

CONTINUE 


K = 0 

DO 85 K = 1 
R(K,K) = 
RINV(K,K) 
CONTINUE 


DO 88 I = 1,2*MODE 
DO 89 J = 1,2*MODE 
FCI) = PHICL;J) 
CONTINUE 
CONTINUE 


Fike PREPARE G MATRIX FOR RICDSD 


CALL MATRAN(B,188,2*MODE ,3,BT,3) 


CALL MATMUL(B,188,2*MODE,3,RINV,3,3,TEMP,NH) 


CALL MATMUL( TEMP, NH, 2*MODE ,3,BT,3,2*MODE,G,NG) 
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Fe Fete te We Fe Fee Fe Fe Fe WK Me FW Fe MK We FCN FW Hee FCW FETE FE FC HE FCC HCFC TE Fe IC He Ie ICH IC WPCA TCH Me He He CT Fe FN I 


BEGIN RICCATI GAIN CALCULATIONS 
FRIAR IAAI AIAN HI HIRI I RICHIE TK IK TKK KICK IK IK IKIKIN 


Glace C2" C) Ca 


CALL RICDSD(NF ,NG,NH,NZ,2*MODE,4*MODE,F,G,H,Z,W,ER,EI,WORK, 


+ SCALE , ITYPE ,IPVS) 
C 
WRITE (6,*) ‘ ' 
WRITE(41,*) ° ' 
WRITE (41,110) 


110 FORMAT (/' THE RICCATI SOLUTION IS:'/) 
DO 120 I = 1,2*MODE 
WRITE (41,1050) (H(I,J),J=1,2*MODE) 
WRITE(41,*) ' ' 
120 CONTINUE 
WRITE (6,130) 
WRITE (41,130) 
130 FORMAT (/' THE CLOSED LOOP EIGENVALUES ARE: '/) 
DO 140 I = 1,2*MODE 
WRITE (6,*) ER(I),EI(I) 
WRITE (41,*) ER(1),EI(I) 
140 CONTINUE 
WRITE (6,150) WORK(1) 


WRITE (41,150) WORK(1) 
150 FORMAT (/' CONDITION ESTIMATE IS: ',D26. 18) 
C 
C dokvetdieker seks COMPUTE GAIN MATRIX L Heke 
C 
CALL MATMUL(BT, 3,3, 2*MODE,H,NH, 2*MODE,TEMP1,3) 
C 
CALL MATMUL(TEMP1,3,3,2*MODE,B,188,3,RR,3) 
C 
DO 103 I = 1,3 
DO 104 J=1,3 
RR(I,J)=R(1,J)+RR(1,J) 
104 CONTINUE 
103 CONTINUE 
C 
CALL DLINDS(3,RR,3,RRINV, 3) 
C 
CALL MATMUL(RRINV,3,3,3,TEMP1,3,2*MODE, BT, 3) 
C 
CALL MATMUL(BT,3,3,2*MODE,PHI,188,2*MODE,L, 3) 
C 
WRITE(41,*) ' ' 
WRITE(41,*) ' GAIN MATRIX L' 
WRITE(41,*) ' ' 
WRITE(41,*) ' ROW 1 ROW 2 ROW 3' 


DO 9155 I=1,2*MODE 
WRITE( 41,1040) (L(J,1),J=1,3) 
9155 CONTINUE 


WRITE(41,*) ' ' 
C 
G PIT Ke Fee Ve TW ese Fe Fe Fee Fe Fe Fe Fe Fe Fe FC FC TC ICICI HCFC FC ICI IOI TCT TPIT FETC FE FETC FONE FTE FACE He Fe FON 
e seeder COMPUTATION OF TORQUES AND COSTS teictoink 
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9000 


gqaqaqnqaqn 


qyvigdagaa, 


210 


MQQAaIIaAIQ QI aa 


FPR HATCH RIC HR FRI IK TKK ICK RICTRICR EK RIKI RTKIKAEERIRRIR IK 


COUNT = 0 
TOTCST = 0. ODO 
TIME = 0.0 


FITCH ICH RII CIERRA RIKER REREREREEEREEEERERREREKREEREEKRIEIKE 


wie SETS LOOP FOR THE NUMBER OF ITERATIONS NECESSARY ***** 


wei TO OBSERVE THE SYSTEM FOR DESIRED LENGTH OF TIME ***%% 
(riniddicidiniitinncnniciiiccnoinickioiccikvivkkctinKqcsvqKrebkkkeckicvitk kirk 


LOOP = INT((MIN*60. 0) /SAMPT) 
PRNT = INT((CMIN*60. 0) /SAMPT) /100. 0) 
PRNTG = INT(((MIN*60. 0)/SAMPT) /2000. 0) 


DO 200 N = 0, LOOP 
TIME = DBLE(N)*SAMPT 


IFCN. EQ. 0) THEN 


IMPLSX = IMPX 

IMPLSY = IMPY 

IMPLSZ = IMPZ 
ELSE : 

IMPLSX = 0. ODO 

IMPLSY = 0. ODO 

IMPLSZ = 0. 0DO 
ENDIF 


Pete Fedede te Fede Tete Fete KF FH TT TIT PTC TTI THICKET KERRIER ERE KERRIER 


wick wikekek CONTROL TORQUE EQUATIONS deactedeteetek tk 
FARIA IAA AA TAT AIHA ARIAT NICHI HI TT THIER ICH IAAI IRIN IR RIN I 


SUM1 = 0. 0DO 
SUM2 = 0. 0DO 
SUM3 = 0. ODO 


DO 210 CT = 1, MODE 
CTADJ = CT + (SMODE - 1) 
SUM1 = SUM1 + L(1,2*CT-1)*X1(CTADJ) + L(1,2*CT)*X2(CTADJ) 
SUM2 = SUM2 + L(2,2*CT-1)*X1(CTADJ) + L(2,2*CT)*X2(CTADJ) 
SUM3 = SUM3 + L(3,2*CT-1)*X1(CTADJ) + L(3,2*CT)*X2(CTADJ) 
CONTINUE 
TCX = SUM1*(-1. ODO) 
TCY = SUM2*(-1. ODO) 
TCZ = SUM3*(-1. ODO) 


week RIK UNCONTROLLED SYSTEM TORQUES KIKKEKKAKEKE 


TCX 
TCY 
TCZ 


0. ODO 
0. ODO 
0. ODO 


FATIH RICK TAKE REIKI ERERRERRRERERIKRERE 


IF(N. EQ. 0) THEN 
WRITE (32,*) ‘IMPULSE X AXIS, IMPULSE Y AXIS, IMPULSE Z AXIS' 
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>) 


oqoaqaiac 


230 


gqaraqan 


WRITE (32,*) ' ' 
WRITE (32,1040) IMPLSX, IMPLSY, IMPLSZ 
WRITE (32,*) 9 * 
WRITE (32,*) 'CONTROL TORQUES TCX, TCY, TCZ' 
WRITE (32,*) ° ' 

ENDIF 


IF (N. LE. 20) THEN 
WRITE (32,2000) TIME, TCX, TCY, TCZ 
ENDIF 


IF (MOD(N,PRNTG).EQ.9) THEN 
WRITE( 30,1036) TIME,X1(7),X1(10) ,X1(30) ,X1(50) ,X1(80) ,X1(.100) 
WRIDE( 3071030) —liME,X107),X1(10),X1( 20) ,X1030), X1(40) ,X1(50) 
ENDIF 


HICK I TTI HIKARI EERE REECE 


KERRI SYSTEM COST FUNCTION CALCULATION KEKRK KA 
HAHAH HA HAI IKI HICK IKI III INIT RIKI RR IRINA RIK I 


SUMC = 0. 0D0 
ENERGY = 0. ODO 
CNTCST = 0. ODO 
COST = 0. 0DO 


DO 230 CF = 7,MODAL 
MODEN(CF) = MODEN(CF)+(X1(CF)**2) *LAMA(CF)+X2(CF)**2 
SUMC = SUMC+(X1(CF)**2)*LAMA(CF)+X2(CF)**2 


CONTINUE 
ENERGY = SUMC 
CNICST = (TCX**2)*RM+( TCY** 2) *RM+( TCZ** 2 )*RM 


COST = ENERGY -+ CNICST 
TOTCST = TOTCST + COST 


IF (MOD(N,PRNT).EQ.0) THEN 

COUNT = COUNT+1 

WRITE(33,2000) TIME,ENERGY , CNTCST, COST 
ENDIF 


Fesete ewe tere Fes Fe se Fe Fete re Fete te te Xe Fe HK Ke IH I HIT HK HK AK KKK KK KKK KK KK 


ieitieiciciei STATE UPDATE EQUATIONS deine 
Fete ke tedeedkesesdedetke kisi ke Keke keke Ke seke Fete Fe TK TK IT ICT RICHI KIKI FF Te ick Ket ke 


DO 220 KA = 7,MODAL 


K = KA-6 
X1T=PHII(1,1,KA)*X1(KA)+PHII(1,2,KA)*X2(KA)+B((2*K-1) ,1)*TCX+ 
ie B((2*K-1) ,2)*TCY+B(( 2*K-1) ,3)*TCZ+BN((2*K-1) ,1)*IMPLSX+ 
ie BN((2*K-1) ,2)*IMPLSY+BN( (2*K-1) ,3)*IMPLSZ 
X2T=PHII(2,1,KA)*X1(KA)+PHII(2,2,KA)*X2(KA)+B(2*K, 1)*TCX+ 
+ B( 2*K, 2)*TCY+B(2*K ,3)*TCZ+BN( 2*K, 1)* IMPLSX+BN( 2*K,2)* 
ft IMPLSY+BN( 2*K,3)*IMPLSZ 
X1(KA) = X1T 


off 


X2(KA) = X2T 


220 CONTINUE 
C 

200 CONTINUE 

C 


WRITE (35,3000) COUNT 
RTOTAL = TOTCST 
WRITE (34,3002) MODE,RTOTAL 


DO 235 K = 7,MODAL 
RMODEN(K) = MODEN(K) 
WRITE (31,3002) K, RMODEN(K) 
235 CONTINUE 


666 CONTINUE 


VAT KIKAIHHKRAR IK RIKI KR RRE 


C 
C 
C Wick CHANGE THE VALUE OF R FOR NEXT RUN??? Wickiek 
Cc KIKAKKEAKIKKIK RIK IIR IE SAKE IAIRAEINETNIK IK IARI IAAP IAI 
C 


WRITE (6,1008) 
540 WRITE (6,*) 'DO YOU WANT A NEW R VALUE? (Y/N) °CAPS ONLY ' 
READ(*,1010) RAGAIN 
IF(RAGAIN. EQ. 'Y') THEN 
WRITE (6,*) 'ENTER NEW R_ VALUE: : 
READ (6,*) RM 


GOTO 550 
ELSEIF(RAGAIN. EQ. 'N') THEN 
GOTO 530 
ELSE 
WRITE(6,*) ‘YOU MUST CHOOSE "'Y" OR ''N'', SELECT AGAIN. ' 
GOTO 540 
ENDIF 
C 
C Jndcccdddicidindck divin inidcicicicinividicticiciicinicciniciccinicincniccikicck 
C Ferrer RUN PROGRAM AGAIN OR QUIT??? HHA 
C Hie ded sete HCH AH eI I He Te HH IIe eve eT TCH RI ICTR AIH HR RII I IR IIE 
C 
530 WRITE(6,*) ‘DO YOU WANT ANOTHER RUN? (Y/N) °CAPS ONLY ' 


READ(*,1010) AGAIN 

IF(AGAIN. EQ. 'Y') THEN 
GOTO 500 

ELSEIF(AGAIN. EQ. 'N') THEN 
GOTO 599 

ELSE 
WRITE(6,*) 'YOU MUST CHOOSE “Y'' OR "N", SELECT AGAIN. ' 
GOTO 530 

ENDIF — 


FEAT CA TACT LET ICHC TCT ATTIC PEPTIC RICH DC ICR ITA CACHAN R AA ARIA IAI ERI 


C 
C 
C FICK FORMAT STATEMENTS bala Ray 
C tidciidddicticinininicicticicceccinincicciciinccicicicincccicicictcicicickiticticttk 
C 


700 FORMAT (' ‘','STARTING MODE NUMBER: ‘',I2) 


701 FORMAT (' ','NUMBER OF MODES SCANNED: ',I2) 
702 FORMAT (' ','NOISE INPUT NODE: ',I3) 
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703 
704 
705 
712 
713 
1001 
1002 
1004 
1005 
1008 
1010 
1035 
1036 
1040 
1050 


670 


660 
650 


oamda a 


FORMAT (' ','INITIAL R VALUE: ',E12. 4) 
FORMAT (' ','SAMPLING TIME: ',E12.4) 
FORMAT (' ','DAMPING FACTOR: ',E12.4) 
FORMAT (' ','LAST CONTROLLED MODE: ‘',I2) 
FORMAT (' ','OBSERVATION TIME: ',F5.1,' MINUTES' ) 
FORMAT( 1X, A6) 

FORMAT( 1X,8E15. 8) 

FORMAT( 1X,//) 

FORMAT( 1X ,60X,E11. 5) 

FORMAT(1X,////) 

FORMAT(A1) 

HORMAG(: ,F7.2,2X,5¢bi2.6,2X)) 

MORMADG ,F7.25uXeoCE11.5,1X)) 

FORMAT(' ',3(E15.8,5X)) 

FORMAT(' ' ,4(E12.6,2X)) 

FORMAT( 1X ,F7. 2, 3X,3(E15.8,3X)) 

FORMAT(' ',T5,E15. 8) 

FORMAT( 14) 

FORMAT(F7. 2,2X,E12. 5) 

FORMAT(1I3,2X,E12. 5) 


STOP 
END 


PePe dee Xe TEPER T CICA RT CIC TCR TET TCT TOR ICTR IERIE TE FOR TOT TE FOPCTET TOR FOP RARER RRR RRR FC 


SUBROUTINE TO MATRIX MULTIPLY 
PAAR IN III I IRI IRIE RTI IIR IK IITA IRIE 


SUBROUTINE MATMUL(M1,LD1,R1,C1,M2,LD2,C2,MP,LD3) 
INTEGER R1,C1,C2,LD1,LD2, LD3 
REAL*8 M1(LD1,1),M2(LD2,1),MP(LD3,1),SUM 


DO 650 I = 1,R1 
DO 660 J = 1,C2 
SUM = 0. 0DO 
DO 670 K = 1,C1 
SUM = SUM+M1(1I,K)*M2(K,J) 
CONTINUE 
MP(I,J) = SUM 
CONTINUE 
CONTINUE 
RETURN 
END 


FOFETA TE TAT CACHE PEPCK CTE R TCPCT ALC FOR POPPER ICV A CHIC PCT ASCH ICSCR AIC ICTCR ACRE ARRTE 


SUBROUTINE TO TRANSPOSE A MATRIX 
FHP A IARI IIE IIR IRIE RIT IR RIKI AA AAA AAAI AAAI AIA AAI 


SUBROUTINE MATRAN(MX,LDX,R1,C1,MT,LDT) 
IRTEGER R1,C1,1,J,LDX,LDT 
REAL*8 MX(LDX,1),MTCLDT, 1) 


DO 680 I = 1,R1 
DO 690 J = 1,C1 


50 


MT(J,I) = MX(I,J) 
690 CONTINUE 
680 CONTINUE 
RETURN 
END 
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APPENDIX B. MODAL DATA 


A. MODAL SLOPES FOR NODE 69 


The modal slopes for the control node are: 


Mode 


OOmn~IA UN & W Pe 


‘oo @4fmeoeornfst#fts#fre°d’ 1 Oe ao rt ' © OO 8 8 4 


, 
X 69 


x 


»142743E-15 
. 148348E-14 
960969 2b 
. 203930E-04 
. 138406E-06 


178740E-06 


. 134835E-04 
. 153043E-05 
. 814600E-07 
soos 550-06 
»l31939E-05 


446746E-08 


76372768 707 
. 485695E-07 
70209505 -0> 
. 164529E-05 
w2623 75-06 
8155525 -06 
so 7O9CsE-07 
"98 22928-06 
- 267182E-05 
1 229879E-06 
SPLIT G MNS 
.412678E-06 
.644907E-05 


-. 189005E-04 


io 1-°Oo°e tO 1°50 ' 4 oO 1 OO ts 


Resp aN 8 oO) 
. 225840E-04 
. VOSE2ZZE=05 
- 263199E-04 
. 648130E-04 
. 228027E-06 
.112004E-04 
. 102071E-04 
- 2/0660E-04 
- 460758E-04 
. 902668E-05 
SOSOcOIE 05 
- 526097E-06 
-574919E-05 
. 202286E-04 
.434762E-05 


rr OS0251h=05 


2) 


» 176625E-04 


rooor 41 © 


roo 1jononcrijinonnnonrs#Fr”nedteos:ie. 


Ooo 1 


Ceo Oo Ot tt OS 414 


x’ 
69, 


.710767E-14 
2127912E=15 
- 174564E-13 
. 264468E-14 
.311234E-04 
- 127699E-06 
-S50029E-06 
S29 305-05 
.578348E-05 
, 51O083Z2E<-05 
. 101384E-06 
. 544769E-08 
2/7 LIPS 0y, 
-424607E-06 
. 266580E-06 
, 160255n-05 
- 300355E-05 
- 198735E-04 
.102157E-04 
» 303975E-04 
.312356E-04 
- 223818E-04 
; L81S74h 04 
706265 n 505 
. 360769E-04 
-493369E-05 
.619046E-05 
72219538 -05 
- 188184E-04 
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B. MODAL SLOPES FOR NODE 23 


The modal slopes for the shuttle docking point are: 
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C. MODAL SLOPES FOR NODE 55 


The modal slopes for the alpha-joint are: 
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